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Section 1 

INTRODUCTION 

This document represents the culmination of a theoretical investigation 
into the geometry changes taking place during the burning of a solid rocket 
motor — especially m the regions of atypical burning surfaces. The final 
destination of these results is intended to be the Boeing developed SRB-II 
solid rocket performance prediction model , and is thus a method of freeing 
that model from excessive geometric input and a reliance on tables of burn 
area as a function of burn distance. It is, however, a general geometric 
codification and can be of valuable use m other models as well. 

The performance of a solid rocket motor (SRM) with a given propellant, 
that is, its internal pressure-time history, is largely dependent: upon its 
internal geometry. Since the propellant itself is a solid, several basic 
shapes can be selected for the interior port. The naming for each configura- 
tion is based on the geometrical figure inscribed by the port on a plane 
taken perpendicular to the motor (longitudinal) axis. Figure 1-1 presents 
the basic shapes m common use today, along with the name and abbreviation 
(if one exists) associated with each. 

In response to the particular circumstances surrounding the SRM committed 
to the Space Shuttle program, it was decided that a slotted cone and cylmdri- 
cally perforated (CP) grain combination would suffice. The internal pressure in 
the SRM is directly dependent on the total amount of propellant surface area 
exposed to active burning. The star (here, the terms slotted cone and star 
are used interchangeably) configuration at the head end (Figure 1-2) will 
allow a large surface area to be eiqposed at ignition, producing a rapid build- 
up of internal pressure and thus a high thrust. Since the propellant there 
will burn so rapidly, it is the general case that this particular region is 
consumed long before the rest of the propellant in the motor. The CP section 
then takes over as the prime gas generating region, not too long after ignition. 
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The analytical model described here was specifically constructed for the 
slotted cone and cylmdrically perforated combination of cross sectional 
geometry. The number of slots m the slotted cone remains variable. Modeling 
a standard star is based only on proper variable assignments and is handled 
in the same way as the slotted cone. For this reason, as indicated above, the 
two will be called "star", because there is no analytic difference. 


Typical surface area traces for a slotted cone region and a CP region are 
given in Figure 1-3. 




SURFACE/OHME TRACE FOR A SURFACE/TIME TRACE FOR 

STAR REGION A CP REGION 

Figure 1-3. SURFACE AREA VERSUS TIME TRACES IN DIFFERENT MOTOR SEGMENTS 


Notice that the star segment is regressive, that is, the exposed surface 
area decreases rapidly with time. Hence, its influence during the burn 
diminishes. Just the opposite is true of the CP region. .Its burning surface 
increases during the burn, and its influence becomes great. For this reason, 
the burning m a CP region is termed progressive. In passing, it is pointed 
out that the pressure trace is directly related to a summation of the two 
traces in Figure 1-3. 

Large SRM's, such as those on the Space Shuttle, require casting the 
propellant in several sections, called casting segments; and these, in turn, 
are joined together (usually at the launch site) to form the motor. The gaps 
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m the propellant caused by the joining of these casting segments are called 
slots. These slots present new surfaces to the burning. 

In order to "program" the burning of a solid rocket motor to conform to 
a given thrust-time history, these slots may or may not have to be inhibited 
to some degree to eliminate or control these new burning surfaces . A non- 
burning ablative material can be placed against the propellant bounding a 
slot to control these otherwise bothersome areas. If no burning is wanted m 
this region, the propellant is covered completely with the inhibitor. If some 
burning is needed, the inhibitor can be built up to leave an annulus-shaped 
region bare to the burning (Figure 1-4) . In such a region, burning not only 
takes place radially, but also longitudinally. 



INHIBITOR 


PROPELLANT 


Figure 1-4. INITIAL PROPELLANT FACE IN A PARTIALLY INHIBITED SLOT 


In most regions of an SRM, burning takes place in the .radially outward 
direction. Special regions do exist m some motors where the burning takes 
place not only radially outward, but possibly m other directions as well. 

The just described slot burning is an example. These are the regions in 
which special care must be taken with any model. The two- and sometimes three- 
dimensional regression of the grain leads to a complex geometrical configura- 
tion. It is for these configurations that this geometric model was developed. 
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The following discussion progresses along well defined lines. Each 
region of a motor (exclusive of the head dome region and the aft dome region) 
will be handled m turn in its own paragraph. The general progress is from 
the simple to the difficult, in hopes that both geometrical intuition and 
analytical sophistication can be fostered. 

The first three sections of the analytical model deal, respectively, with 
the CP, star, and partially inhibited slot regions. The final, and by far the 
most involved, section describes the region of transition between the star and 
CP segments, where burning proceeds m three dimensions. A section of verbal 
description and general mathematics precedes the analytic model , and more 
mathematics of a less influential nature is presented m two appendixes . Con- 
tinuity can be maintained by passing over the mathematics section if the 
equations in later sections are to be accepted without argument. 

Some results for the Shuttle- committed SRM gram design are given m 
Appendix C. 
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Section !i 

SOLID ROCKET GEOMETRY 

This section gives a verbal description of the burning of the propellant 
in an SRM. Also included is an extensive description dealing specifically 
with the mathematical methods which will be used to model the gram regression. 
The notions of increment dividing planes and reference planes are introduced, 

„ and the philosophical differences between the two are explored. 

Care should be taken to absorb paragraph 2.4 which deals with the volume 
and surface area formulas used within the model. Those formulas, along with 
many others within the body of this work, have been developed from information 
taken from References 1 and 2 (see, however. Appendix A) . 

2.1 ATYPICAL GEOMETRY AND GRAIN REGRESSION 

Any portion of the gram which is not interior to a CP or star segment is 
likely to develop an atypical (complex) geometry during the burn. These areas 
are characterized by the fact that the geometry makes a radical change; for 
example, a portion of the CP becoming star shaped. In this sense, regions of 
complex geometry occur in the gram surrounding a partially inhibited slot and 
in the region of transition between star and CP, The portion of the motor 
which is uniformly star shaped does not, for example, constitute a region of 
complex geometry, since a six pointed star remains that until the burning has 
consumed the entire region. 

Even though the star section is not, per se, a region of atypical geometry 
the region has been mathematically modeled here in order to be compatible with 
the three-dimensional burning taking place in the star to CP transition region. 
This also requires a new modeling of the burning in the CP region. In so doing 
the entire surface of the grain can be handled with analytical methods, through 
functions depending only on total distance burned. 

Figure 2-1 shows the geometry changes which can take place m different 
regions of the motor. Together with those shown, the transition between star 



NORTHROP SERVICES, INC. 


TR-1726 


\ 


/ 


STAR 

SEGMENT 


/ 



\ 



CP - SEGMENT 


CROSS SECTION TAKEN PERPENDICULAR TO CENTER LINE 



CROSS SECTION TAKEN PARALLEL TO CENTER LINE 

Figure 2-1. GRAIN REGRESSIONS IN DIFFERENT PORT CONFIGURATIONS 
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and CP segments offers a challenging exerco.se in planar renditions of spacial 
regressions. Some three-dimensional offerings are given later in the theore- 
tical modeling of that region (Section V) . Three basic assumptions govern 
the burning and produce renditions like those in Figure 2-1 : 

1 . Propellant burning during ignition and steady state occurs along the 
outward normal at each point of the port-gram interface . 

2. Initial arcs of propellant reaching into the port progress to points. 

3. Outward angles progress to arcs. 

The first assumption provides the vehicle for the analytical description, while 
the second two provide means by which analytical discontinuities , whether real 
or potential, may be handled. 

Pencil and paper gram regressions are made by taking the original cross- 
sectional configuration and drawing curves toward the case wall at equal 
distances from the initial curve. It is not unreasonable, then, to have 
assumptions like 2 and 3 above. 

Before turning to the analytical description itself, it will be of interest 
m later theoretics to go through some general, but useful, developments. The 
next section is devoted to this. 

2.2 GENERAL MATHEMATICS 

It will be of interest further on to have a few general notions adequately 
developed. The information contained m this subsection will be called when 
needed by appropriate names posted m this paragraph. This paragraph may be 
passed over without disruption of reading continuity, except for the definitions 
herein presented. 

Two mathematical notions share equally in importance in this geometry 
model. The first is in response to seeking what happens when families of 
concentric circles intersect families of parallel lines, and the second is a 
method for analytically describing familiar shapes m space. 

Consider first the situation appearing m Figure 2-2. On the one hand, 
there is a circle centered at point (a, b) in the x-y Cartesian coordinate 
system, and on the other, a straight line, whose equation is given. 

2-3 REPRODUCIBILITY OF TH& 

original page IS POOft 



tvflJi 





Figure 2-2. DEFINITION OF A CIRCLE-LINE PAIR 


Assume that the radius of the circle grows at a specified rate, while the 
original line moves in the direction of the circle center at the same rate . 

A circle-line pair is defined to be the circle centered at (a, b) with 
radius r + Ar along with the line which has moved Ar distance from the original 
line, m the direction indicated. In this respect, the circle and the line 
°f Figure 2-2 represent a circle-line pair. The distance Ar is taken to be 
any real number such that Ar > -r. 

Considering only the points of intersection (if extant) of circle-lme 
pairs for all Ar >_ -r, what is the functional analog of this curve? (Clearly a 
curve results because of the continuous nature of the movement.) 

The original circle has equation 
(x - a) 2 + (y - b) 2 = r 2 
and the original line has equation 
ax + By + y = 0. 
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Suppose, at one point, the radius of the circle has grown to r + Ar. The new 
circle has equation 

(x - a) 2 + (y - b) 2 = (r Ar) 2 


If <b represents the angle of inclination of the line with respect to the 
x-axis {that is, <}> = arctan (a/B) ) , then the new line has equation 

ax + By + y + Ar/cosd> = 0 


Thinking of the above two equations as equations m three variables (x, 
y, Ar) , one variable m the circle equation can be eliminated using the line 
equation. For example, eliminate Ar. Solving the line for this quantity, 
obtain 


Ar = -{ax + By + Y) cos6 
or 

Ar = - (ax + By + Y) 


where a = acos<f>, and so on. Substitute this value into the circle to obtain 


(x - a) 2 + (y - b) 2 = (r - (ax + By + Y) ) 2 . 

Expanding and rearranging the form 

(1 - a 2 ) x 2 - 2aB xy + ( 1 - Q 2 ) y 2 - 2 (a - ra + ay) x 


- 2 (b ~ rB + By) y + a 2 + b 2 - r 2 ° — - - 2 


2ry - y =0 


is obtained. If these assignments are made: 


A = 

i ~2 

1 - a 

D = -2 (a - ra + ay) 

B = 

-2aB 

E = -2(b - rB + By) 

C = 

1- B 2 

2 2 2 
F = a + b - r + 


( 2 - 1 ) 


the equation has the form 
2 2 

Ax + Bxy +Cy + Dx + Ey + F = 0 

Notice that this is the general second degree equation. It then represents 

♦ “ 2 
either a hyperbola, a parabola, or an ellipse, depending on B - 4AC > 0, = 0, 

or < 0, respectively. But 
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B 2 - 4AC = 4 a 2 6 2 - 4(1 - a 2 ) (1'- g 2 ) 
or 

B 2 - 4AC = 4 (a 2 + I 2 - 1) 

2 

The conditions on B - 4AC are then translated into conditions on the line 
coefficients a and g. These results are summarized as: 

Hyperbola: 

Parabola : 

Ellipse: 


The family of points arising from the intersection of circle-line pairs 
represents a general second order polynomial m two variables. Specifically, 
this second order equation is a hyperbola, parabola, or ellipse depending only 
on the location of the pair (a, g) with respect to the unit circle — "outside," 
"on," or "inside" can be translated as "hyperbdla, " "parabola," or "ellipse," 
respectively. 

These particular shapes play important supportive roles during the model- 
ing of a burning SRM. It will be remembered that the initial grain surface, 
for the most part, is piecewise linear. At abrupt angular changes, burning 

I 

causes arcs to emerge, which grow at the burning rate. These circular arcs, 
then, are adjacent to linear regions of gram surface; and during burning, the 
point of intersection of the grain surface with the formed circular arc causes 
one of these second-degree equations to become "law of motion," used to track 
that point of intersection. Predictions can then be made as to the location of 
a circle-line pair intersection at important times, such as end of bum. For 
that reason, such second-degree formulas will be termed "burn -back" curves — a 
bum-back parabola, burn-back hyperbola, or burn-back ellipse. 

Now, turn to the other major notion — that of finding an analytical formula 
for familiar geometric objects in three space. Here, a powerful method of 


B - 4AC > 0 


—2 —2 
a + g > 1 


B - 4AC = 0 


a 2 + g 2 = 1 


B - 4AC < 0 


~2 —2 
a + g < 1 
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advanced calculus* is used to obtain specific equations. The method ultimately 
rests only on the fact that in space, there are two planes which intersect at 
right angles . 

Two examples will serve to illustrate the method. The first develops the 
equation of a specific torus m space; the other develops the' equation of a 
cylinder in space which leans into the first octant with specified angles of 
inclination . 

EXAMPLE 1: TORUS 

In an x-y -2 Cartesian coordinate system (Figure 2-3) , assume that a circle 

centered at (a,b,0) of radius r is revolved about the line x = c. The circle 
2 2 2 

(x - a) + (y - b) = r is called the base curve; the other is called tne gen- 
erator curve. The generator is a circle centered at (c,b,0) of radius a - c. 
Revolve the circle as described to obtain a shape like that in Figure 2-4. 

The base circle can be expressed as a multiple valued function of y as 
x = f(y) = + yj r 2 - (v-b) 2 + a 

Choose a point P with coordinates (x = f(y ), y , 0) on the base circle, 
o o o o 

Choose an arbitrary point P, (x,y,z) on the circle that P makes when revolved 

1 o 

about x = c. Since the revolution takes place parallel to the x-z plane, then 
y o = y^; x^, and z satisfy the circle equation 

2 2 2 
( X . - c r + z: = (x - cr 

1 1 o 

Substituting appropriate quantities, obtain: 

(x - c) 2 + z 2 = (f(y ) - c) 2 
1 1 o 

(x 1 - c) 2 + z 2 = (f(y x ) - c) 2 

, ,2 L 2 , , / 2 ' , . 2 ,2 

(x.^ - c) + z = ( ± yj r - (y 1 - b) + a - c) 

*See Reference 2 3 Pages SI 5-6 20. 
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There is nothing magical about the subscripts, and they can be dropped to 

obtain 

(x - c) ^ + z 2 = (± y r 2 - (y - b) 2 + a - c) 2 

which is the appropriate formula for the torus . Choice of the + or - sign 
separately yields the upper and lower regions, respectively. In summary, the 
equation of a torus with generator curve, circle centered at (c,b,0) radius 
a - c and base curve: circle centered at (a,b,0) -with radius r, is given by 

2 2 1 2 2 2 

(x - c) + z = (± y/r - (y - b) + a - c) 

EXAMPLE 2: AN OFFSET CYLINDER 

Consider a cylinder as drawn in Figure 2-5, whose centerline rests m the 
x-z plane. 



Figure 2-5. GENERATION OF AN ASKEW CYLINDER 
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This cylinder will be described first, then the figure will be rotated into 
the first octant to obtain final resultant offset. 


Using the same methods as in the last example, the base curve is a 
straight line indicated by &, while the generator curve (the generator curve 
must be in the x-y plane in Figure 2-5) is an ellipse given by 

(x - h) 2 . (y - k) 2 

+ — = — - — = 1 

2 ,2 

a b 

and has the following attributes : 
a = r/sina 
b - r 


e = cosa, e - eccentricity 

Next, choose a point P q (x q , Y 0 ' 0) on base curve. The base curve has 
equation 


x = x + t cosa 
o 


y = = 


■ = b 1 -/x - h\ 2 + 

° y (-V-) 


..., all real numbers, t. 


z = t sina 

Since x and z are the only two equations containing t, t can be eliminated as 

x = x + t cosa 
o 


z = t sina 


x - x. 


t - z/sma - 


cosa 


so that 


(x - x ) sina *» z cosa 
o 


or 


But ' 


z = (x - x ) tana 
o 


x o “ 3 - /l - 




+ h, y » y 
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Notice that when z = 0, there is the cross section (trace) 


M-M 


which is the generator line. When a = 90°, then the right circular cylinder 
is obtained, given by 


2 2 2 2 
(y - k) + (x - h) = a' = b 


Taking a plane parallel to the x-y plane at z = c, there is 

(^) 2 * g ) 2 - 1 

which is an ellipse centered at (h + ccota, k, c) with eccentricity as 
described. 


Hence, the above formula is truly that of the described cylinder. Next, 
twist the cylinder m the positive y-direction to move it into the first 
octant (Figure 2-6) . 



Figure 2-6. AM OFFSET CYLINDER 
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It is easy to see that equation (2-1) must be updated to the form 
/ y - (k + zcosBcota) V + / x - (h + zsmgcotg) ^ 


(2-2) 


Notice if a = 90 ° , there is 

,2 ,\ 2 


M * (Hr 1 ) ■ 


= 1 a = b 


If 0 = 90' 


(y - k\ 2 . /x - (h + zcotcO^ 2 

nr/ + (- — ^ 1 = 1 


and all things reduce to the last case. The trace m the plane z — c yields 
/y - (k 4- ccos8cota )\ 2 / x - (h + csm8cotct )\ 2 _ 

\ b / * V a / 


which is an ellipse centered at 

I 

(h + csingcota, k + ccosgcota, c\ 
namely, a rigid translation of the generator curve to that point. 


In summary, it can be stated that the equation of a cylinder skewed from 
the y-z plane by angle 8, and from the x-y plane by angle a, and having radius 
r must have equation 

2 

( y - (k + zcosBcota n + / x - (h + zsingcotcQ j _ ^ 

b ; \ * / 

where 

b = r a = r/sina. 


For an application of both the cylinder and the torus equation, see the 
section on modeling the transition region. Section V. 

2.3 REFERENCE PLANES AND INCREMENT DIVIDING PLANES 

Initial geometry is conveyed to the model by way of a number of reference 
planes . These are simply planes positioned perpendicular to the centerline of 
the motor at user-defined intervals. They convey such information as the 
position along the motor axis from a given reference, initial port radius. 
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and other data to specify the distinctions between star and CP. These planes 
are also used to convey insulation thicknesses by specifying the radius from 
the centerline to the case wall. 

Since only a few reference planes are needed to adequately define an 
internal geometry, it is then the general case m a large motor that these 
planes are wide apart. An internal ballistics solution using only the 
reference planes, then, cannot have needed accuracy, so a shrinkage of the 
world view is necessary. This is accomplished by the device of the increment 
dividing plane (IDP) . 

In order to space geometric data closer rogether than is possible using 
a small number of reference planes, the increment dividing planes are intro- 
duced. Since these are set at positions in the motor by the code, an other- 
wise lengthy and necessarily redundant list of reference plane geometry input 
constants is reduced to a bare and necessary minimum. 

The increment dividing plane becomes the ma^or vehicle to transmit 
geometric data to the code. The values at any increment dividing plane include 
those at reference planes with the addition that on this level, total burn 
distance is tracked. The increment dividing plane is not a static notion, as 
is the reference plane, and some increment dividing planes actually change 
position to track radical changes in geometry. 

Information is conveyed to the increment dividing plane level by first 
positioning each IDP at reasonable intervals between reference planes, making 
sure that each reference plane position is also occupied by an increment 
dividing plane. Other values are transmitted to the increment dividing plane 
using linear interpolation of like variables of the two-reference plane 
bounding a given IDP and the position of the IDP. At this point, attention 
is called to the necessity of adequately placing a sufficient number of 
reference planes along the motor axis to convey important changes in taper and 
in insulation thickness . The interior grain configuration is guaranteed to 
change linearly in the areas mentioned, so the prior level of generality is 
maintained by the use of linear interpolation routines. 
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After the increment dividing planes have been set and geometry input has 
been conveyed to this level, the reference planes are discarded, and only 
serve the code as a repository of information. The dynamic changes in geometry 
are calculated on the increment dividing plane level , and the reference plane 
concept is reduced to one of only transmitting initial data. No calculations 
are ever made on the reference plane level. 

Since the planes are perpendicular to the centerline of the motor, the 
ideas of port fuel perimeter and port area make sense . These are calculated 
and tracked at each increment dividing plane, with this information at two 
adjacent IDP's, using volume and surface area formulas yet to be described, 
the port volume and burning propellant surface area can be calculated. These 
values represent the ultimate reason for the existence of this model, since 
these values are, in turn, fed to the internal ballistics calculations, and 
these govern the regression of the gram. 

2.4 VOLUME AND SURFACE AREA FORMULAS 

This paragraph describes the philosophy and use of the volume and surface 
area formulas used m this geometric code. Before turning to these, however, 
the restriction necessary on the model is restated, since the formulas depend 
on this. 

It is assumed at the onset that the regression of the grain occurs normal 
to the propellant surface at any time. With this assumption, and the fact that 
initial surfaces vary piecewise linearly along the motor axis, it would seem 
that the best formulas are the family of cone frustum mensuration formulas . 

With these formulas, gram linearity between increment dividing planes is 
tacitly assumed. (Notice that linearity between dividing planes is stated 
and not reference planes. The burning itself may cause a nonlinear burning 
surface to emerge, which is approximated by a piecewise linear surface, with 
changes in taper taking place at IPP’s. This will admit a better approximation 
than one between reference planes.) 

However, in certain areas of complex geometry, even frustum formulas are 
not adequate, so a new volume formula is introduced which is particularly 
suited to the analytic nature of this code. 
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Figure 2-7. A CONE FRUSTUM 


Surface area (lateral) 

S = 7T + r 2 ) ^(r~ - r 2 ) 2 + h 2 

Volume 

V ■f.hdj 2 + r 2 2 + rjr 2 ). 

From the geometric routines, port perimeter and port area are determined. 
Consequently, with these known, the surface area and volume equations should 
be compatible with that knowledge. A digression is made here to render these 
formulas so compatible. 

Surface Area 

S = tt (r^ + r 2 ) 

= (irr 1 + -rrr 2 ) 

manipulations will yield 

S = (2irr., + 2irr ) J (2irr - 2irr„)^ + 4ir^h^ . 

4tt 1 2 V 1 2 

But 2-irr^ = P^, 2irr 2 ~ V ^ — the perimeters at top and bottom — (Figure 2-7) so 
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S = h ®1 + P 2 > V (P 1 - * 2 > 2 + <2,rh > 2 • 

Volume 

V - I' 1 * ( ' r x 2 + V + r l r 2> 

1/2 2 / 2 p - * 

= J h ( + Ttr 2 + (m^ ) (Trr 2 ) 

2 2 

If = irr^ , A 2 = trr 2 , the port areas at top and bottom, then the volume 
formula can be written as 



With these equations in proper form, attention is turned for the moment 
to the irregular geometry occurring m the transition from star to CP. The 
grain itself may vary nonlinearly in this section, so that frustum formulas 
may yield only approximations to the true values. This can be circumvented 
m the case of volume in the following way: 

A general prisma toid is defined to be a solid whose area, A , at a 

2 

section parallel to and distant z from a fixed plane can be expressed as a 
polynomial in z of degree not higher than three. If is the area of a lower 
base, B 2 the area at an upper base, and M the area at the midsection, then 

V = ~ h (B 1 + 4M + B 2 ) . 

Examples of general prisma toids include cones, cylinders, spheres, and spheroids. 

This particular formula presents a powerful tool for the calculation of 
the volume in an irregular, nonlinear region. It is a generalization of the 
frustum formula just given, but does not depend on linearity. It does, however, 
need the extra calculations for M, the area of a midsection, and for this reason, 
the general prismatoid formula is used only when necessary. 



Unfortunately, there is no such generalization for the surface area formula. 
Instead, the formula for a cone frustum is used repeatedly where the interval, 
between IDP’s is broken into several (equal) smaller subintervals. This 
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piecewise linear approximation is as close to true values as is computationally 
feasible and expedient. It must be remembered that values of port area and 
port perimeter are available at any position along the motor centerline by 
virtue of the functional nature of the geometry (it -depends only on total burn 
distance) ; and so the breaking up of the interval for surface area or volume 
is only an exercise in function evaluation. 

The formulas for surface area and volume are different still m the 
regions of partially inhibited slots. There, the volumes and surface areas 
are described analytically through use of the calculus and surfaces and volumes 
of revolution. This will be handled in the section on partially inhibited 
slots (Section IV) . 
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Section 1 1 1 

NORMAL GEOMETRY 


Although the homogenous CP and star regions of the motor are not regions 
of complicated geometry, and very efficient codes already exist for their 
modeling, this section offers still another mathematical interpretation of 
the burning in these sections . This was considered expedient to facilitate 
computational continuity at the interfaces between these regions and those of 
atypical geometry. 

Included in this section are paragraphs dealing with specific computa- 
tional quantities m the CP and star regions, m that order. The idea of a 
bum-back parabola is introduced in the star paragraph, and extensive use is 
made of this simple but powerful notion. 

The star algorithm presented here is not m its final and most usable form, 
and will be further updated m the section on the transition region, Section V. 

3.1 CP BURNING 

Consider the view presented in Figure 3-1, that is, a segment of propel- 
lant between two increment-dividing planes, I and 1+1, which is totally within 



Figure 3-1. A SECTION OF PROPELLANT IN A CP SEGMENT 
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the CP region of the motor, and will not change its cross-sectional geometric 
shape. The port volume contained between the two planes is called a mass 
addition region. 

In the general case, the initial port radii at the two positions are 
dissimilar, so that the region between the planes forms a cone frustum. 

Further restricting the view, consider only the cross section at plane I 

(Figure 3-2) . The initial port radius is r, and burn distance is s^. The 

port fuel perimeter and port area at plane I (Figure 3-1) are given at this 

burn distance by 

Perimeter = P_ = 2ir (r + s ) 

I b 

, 2 
Area = A,. = m(r + s, ) . 

I b 



Figure 3-2. AN ANALYTICAL VIEW OF A CP CROSS SECTION 
If the same quantities have been established at plane 1+1, then 

(p i + E W \/ (p i - p i + i > 2 + l2 ' h>: 
Surface Area = SA = 


4tt 
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and 

Volume - V - i h (Al + 

where h is the known distance between the two planes. 

Further refinements to this rather simple model have been made for 
purposes of flexibility and accuracy. The first is to consider the port not 
in an overall view of 2rr radians, but on a sector basis. 


If work is done in this mode, the number of sectors, n, is input to the 

program; and an angle, a, is calculated as 27r/n. Then the perimeter and area 

equations become 

P = a (r + s, ) 

I a 


and the same frustrum formulas can be used to establish the surface area and 
volume of the mass addition region. 


The other addition to the model constitutes a way to treat the final 
burnout in the mass addition region. In this respect, consider the situation 
presented in Figure 3-3. Notice that the burning has caused complete gram 
consumption at plane I; while gram remains at plane 1+1. The propellant 
involved in this region has, of course, a linear surface. In Figure 3-3, 
only the propellant is shown — the case and insulation have been omitted. 


The situation m Figure 3-3 leads to an abstraction like that m Figure 3-4, 
which shows a cross section taken through the centerline perpendicular to 
planes I and 1+1. Figure 3-5 is a further abstraction for the purpose of 
analysis and nomenclature. 

The slope of the grain, M Q , is given by 
(r + s, ) - (r" + s') 


where D is the undirected distance between plane I and 1+1. 
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Figure 3-4. AN ANALYTICAL VIEW OF A CP BURNOUT 
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Figure 3-5. DEFINITIONS OF TERMS USED TO DETERMINE A CP BURNOUT 


The slope of the case, M^, is likewise given by 


M = 
c 


r „ , r 
c ,1+1 - c,I 


An exercise in the use of the laws of cosines yields the position of 
point P relative to the foretangent plane as 

r - r" + s' 

Position of P = Position of I + — — 

M - M 
G c 

At point P, the port perimeter is maximal, and the port area is maximal. The 
volume of the mass addition region is given by calculating two volumes, from 
I to P and from P to 1+1; and the surface area is given by using the frustum 
surface area formula for the region between P and plane 1+1 . Burnout 
occurring from the other direction is handled by the same formula. 


The above two additions to the CP burning routine make that routine a 
general, self-contained code useful m asymmetrical gram studies and without 
need for augmenting routines to calculate mensuration quantities for a partially 
burned-out mass addition region. 
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3.2 STAR BURNING 

Attention is now turned to the star segment of the motor (Figure 1-2) . 
Consider the situation shown in Figure 3-6.. A regression of the grain is 
shown m F-igure 3-7. It is evident (Figure 3-6) from the symmetry involved 
that only one-half of one star peak need be modeled. t The port fuel perimeter 
and the port area for the whole cross section are then integer multiples of 
those quantities involved m the restricted view. 


O 




Figure 3-6. STAR CROSS SECTIONS AND SECTOR DEFINITIONS 


The cross section shown in Figure 3-6 is taken perpendicular to the 
centerline so that the vertex of angle lies on the motor axis. The 

valley side is, in the general case, slanted. Line £ (Figure 3-7) is con- 
structed perpendicular to the valley side through the fillet center represented 
by point Q. 

In Figure 3-7 , the actual calculation of port fuel perimeter and port area 
can take place on an even more restricted view. In particular, four areas of 
distinct geometry present themselves. These areas are: 
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1. The concave portion of the star peak 

2. The linear portion of the valley side 

3. The arced portion of the valley side, emmanatmg from the fillet 
circle . 

4. The linear valley floor. 


O 



The discussion now covers the delineation of these four subregions and 
the calculation of the pertinent quantities m each. This is done by tracking 
a few points along predetermined paths. In effect, the assumption of perpen- 
dicular burning establishes these paths and makes such an analysis rather 
simple . 

The specific points that will be tracked are shown m Figure 3-8 . The 
figure represents two times during the burn, the initial perimeter and one 
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after a distance has been burned. Note that the four points P-j^ through P 4 
are at initial positions, and P^ through V' are at terminal positions. The 
regions for which these points are vertices comprise the sections referred to 
above . 



Figure 3-8. NOMENCLATURE USED TO DEFINE THE BURNING IN A STAR SECTOR 


It is the general case where line P^ ? 2 in Figure 3-8 is not parallel to 
the valley centerline, OC^ . Because of the distances involved in large SRM’s, 
special care must be taken in any analytical routine to use the existing value 
of angle $ p p P. In particular, as seen in the following discussion, the 

X M 

length P' V>' varies trigonometrically with this angle; and at the beginning 
x x — 

of the burn, when a large length of side exists, even this diminutive quantity 
can have substantial effects. 
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It is also seen that the positxon of varies dxrectly wxth this small 
angle. In Figure 3-8, point P xs calculated and serves as an internal 
reference to the program. A small change in the angle mentioned will show 
xtself in the posxtionxng of P. 

Some defxnitxons are appropriate here, before progressing to actual 
theoretical development. The point 0 of Figure 3-8 is on the motor centerline, 
and points 0, C^, determine the plane of the cross section made perpendicu- 
lar to the motor axis. Points and C rest on the case wall. The line OC^ 
is the centerline of the star peak, and the line OC^ is the centerline of the 
star valley. 

Notice that the star peak is concave in nature, the valley sides slope, 
and that at the bottom of the valley, both a small fillet and a linear length 
exist. The valley side is assumed tangent to the fillet as is the line at the 
trough floor. The judicious choice of values for these constraints permits 
the modeling of both a true star or a slotted cone. 

Order can be obtained by imposing a coordinate system, as is done in 
Figure 3-9, where a Cartesian system is selected over a polar system because 
this leads ultimately to a simplification of certain important formulas. One 
of these is the notion of a bum-back parabola'. 

As the bum progresses, notice (Figure 3-7) that the concave portion of 
the star peak proceeds down the star centerline while, at the same rate, the 
valley side progresses along line QP in direction P. If £ represents the 
angle of inclination of the star valley side with the x-axis, then the linear 
equation representing this side can be expressed as 

y = -xtan<t> + (tan^sini + cos$) s. + C 
where represents the distance burned, and 

C = r„ + r_ c os<j> + (r_ + m sin4>) tan<J>) 

4 X Zt 

where r f is the fillet radius m the trough, 

is the length of line on the trough floor, 
and r 2 as indicated in Figure 3-9, is the valley floor radius, 

minus the fillet radius. 
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Figure 3-9. QUANTITY DEFINITIONS IN A STAR SECTOR 


At any time before the annihilation of the concave star peak, the 

equation of this star peak is given by the circle equation as : 

2 2 2 
x + y - (r 1+ s b ) 2 . 

The intersection of the star valley with the concave region of the star 
peak is the point R, (R^ and R"') m Figure 3-7 for various distances 
burned, s b , s b , s~. The means to predict where the concave portion annihi- 
lates is now at our disposal. 

Placing the line equation in a form so that the coefficients of s b in 
both the line and the circle equations, above, are identical, we have 

y 4 xtancft - c _ 
tan<fsin<j> -f cos<£ ~ S b 

Now, in accordance with subsection 2.2 

— _ tan<fr 

a tan$sin$ + cos<J> - sin< ^ 
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H — , . . . , , — y 

tancj>sin<f> + cos 9 

Y = - Ccos<j> 

— 2 — 2 
ct + B = 1, 

and point R of Figure 3-7 tracks out a parabola as it passes through R" and 
R'^ to intercept the star peak centerline. A parabola of this sort — that 
tracks a particular point — will be called a burn-back parabola. The burn-back 
parabola is pictured m Figure 3-10 for the case just worked. 


Figure 3-10. A BURN-BACK PARABOLA IN A STAR SECTOR 

Furthermore, since a, 3/ and y are all known from the line equation, the 

bum-back parabola has the form 
2 2 

Ax + Bxy + cy + Dx + By + F = 0 
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where 


A = 1 - a 2 

c = l - 3 2 

E = -2g (y + r x ) 

B = -2a J 

D = -2a (-Y + r x ) 

— .2 
F = -(y + 


2 

with B -4AC = 0 (guaranteed) . In descriptive form, the parabola is one with 
focus at (0,0) and axis of symmetry the line y = xcot<J>. 

Now, the solution alluded to above — that of a point where the burn-back 
parabola intersects the star peak centerline at E (Figure 3-10) — is the 
simultaneous solution of the parabola equation with the equation of the star 
peak centerline. That equation is given by 
y = xtanifi 

where $ = ir/n, n = number of star points. 

Substituting this expression into the bum-back parabola yields 
2 2 

Ax + Bx (xtanij/) + C(xtan^) + Dx + E (xtanifj) + F = 0 
Manipulations result in the solution m terms of x as 

sm(4> + ijO +_ 1~1 

X = (r - r - r 2 sin^ - r^cosiji) cosijj 

1— cos ($ + >J>) J 

where choice of sign is dictated by the fact that x is to be positive. 

Since the bum-back parabola- is in essence a function of three variables 

(x,y,s, ), the same parabola can be expressed in terms of x and s. instead of 
b b 

x and y. Doing this, and knowing s^, the x-coordinate of point R can be 
evaluated. Hence, there is all the information necessary to know the exact 
position of R during that part of the bum in which the concave portion of the 
star peak is extant. 

In the manner of the above exposition, the bum-back parabola in terms 
of x and s, can be obtained as 

D __ 

Ax 2 + B xs, + Cs, 2 + Dx + Es. + F = 0 B 2 - 4AC = 0 
b b b 

where the letters A-F here have no relationship with those above. In fact, 
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A = sec 


B = - 2sinijisec <j> 


D - - 2 tan<p (r + r^/cosif) + r^tan^) 


E = 2 { (r^ + r f /cos<j> + r^tan^J/costj) - r^) 


C = tan 4> 


F = (r 4 + r f /cos<j> + r 2 tan<jj) - r^2 


Note that indeed B - 4AC = 0 since 

- 4AC = -4sin. 2 gsec% - 4 sec^tan 2 ? and tan 2 $ — sirr ~<psec <j> 


Knowing at any time yields 

-K + D ) * I 


(2BD - 4AE) S, + D - 4AF 
b 


X = 


2A 


where the plus sign was chosen to keep x positive. 


Referencing Figure 3-11, once point R has been placed on the burn-back 
parabola, then the .quantities of port area and port perimeter can be calculated. 
This is done in four segments, as indicated. Knowing the coordinates of R, 
then angle n can be determined. This angle bounds region I. The port perimeter 
contribution from this section is given by 


H (r^s ) 


and the port area m this section as 

1 , ,2 

2 11 (r i + S b> * 


Angle £RSU is known to be ir/n-n, since line US is parallel with line OZ, the 
valley centerline. Consequently, point S is known as the simultaneous solution 
of the set 

= r 4 

| y = xtan <*RSU) 

To derive the contributions from region II, U is taken as a fixed point, and V 
is marked off s^ + distance from U along line UP. Hence, v is known, and 
lengths of line segments RS, SU, UV, and VR are caluclated. The port perimeter 
in this section' is length VR. The Bretschneider-Heron formula for the area of 
the general quadrilateral gives the port area contribution from this section 
(Appendix A) . 
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y - 



Angle £VUW is tt/ 2 - <j>, so both the port perimeter and port area in region 
HI ste known. Length WZ is simply r^ , so the port perimeter and the port area 
contribution m region IV are known, using Bretschneider' s formula. 

The above calculations yield the total value of port perimeter as 

perimeter = p (r + s, ) + VR + (tt/2 - <f>) <r_ + s. ) + r„. 

1 b f b 4 

A similar summation gives the port area. 

After the concave region of the star peak has been completely consumed, 
a situation like that in Figure 3-12 appears. 
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Figure 3-12. A STAR SECTOR AFTER CONCAVE PEAK CONSUMPTION 


Appearing is an analogy to that m Figure 3-11, yet rpO (that is, region 
I has become degenerate) so that there is no need for calculating mensuration 
quantities in that section. The difference does manifest itself, though, m 
the calculation of the positioning of point R. it is now receding down the 
peak centerline, off the burn-back parabola, at an accelerated rate. Hence, 
when R reaches E along the path indicated, R is thereafter positioned as the 
simultaneous solution of the star valley side and the peak centerline . 

The valley side and the peak centerline are represented by the sex 

I y = -xtan* + s /cos* + C 
b 

y = xtarnj 

where C, as before, is r^ + r^tan* + r^/cos*. This indicates that R has 
coordinate representation as : 

R ~ (z — — r (s, /cos* + C) , - — ~^v — t (s. /cos* + C) ) . 

\tan* + tan* b tan* + tan* b / 
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With R known, calculations proceed from region II through IV to obtain 
port perimeter and port area. This procedure is maintained until R reaches 
point P or V reaches point P, as they both reach point P simultaneously. 


While R is progressing from E to P, burn-out m the star valley has 
occurred, so that WZ no longer contributes to the port perimeter and no 
longer is there a circular sector of fr/2 -<j> radians. Point W (Figure 3-13) is 
tracked up the case wall, away from to C^, as the simultaneous solution of 
the two circle equations 


2 , 2 2 
x + y = 

(x - r 2 ) 2 + (y - r 4 ) 2 = (s b + r f ) 2 


The decrement of the sector angle can then be tracked and the perimeter 
measured. 


Figure 3-13. A STAR SECTOR WITH VALLEY FLOOR BURNED OUT 
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In Figure 3-13, angle ^WUV can be measured and, consequently, mensuration 
quantity contributions from area III^ calculated. In area III ^ , a bit more 
work is necessary to obtain the port area contribution (Appendix A) . 

When R passes P, then Figure 3-14 represents what fuel remains in the 
sector. Again, angle £WUR can be calculated, since W is known in the method 
above, and R now is the simultaneous solution of the line and circle: 


y = xtanifi 

(x - r 2 ) 2 + (y - r 4 ) 2 = (s fa + r f ) 2 



V is no longer of interest and is not calculated. Instead, the Bretschneider- 
Heron formula still works for sides RS, SU, and UR to find the port area contri- 
bution from region II. From angle £WUR, the port perimeter m III the only 
fuel left is calculated. The area contribution from lll 2 is updated. 
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These calculations are used to track the rest of the burn — burn-out 
occurring when R reaches 

The above description is valid only in the region of the star and not m 
the star-CP transition region. The above procedure will be augmented when that 
region is assessed. 

Port volume and burning surface area in mass addition regions are calcu- 
lated using the frustum formula described above. 
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Section IV 

PARTIALLY INHIBITED SLOTS 

A partially inhibited slot exhibits irregular geometry during the burn. 
Not only is there a burning component m the radial direction, but also one 
in a longitudinal direction. The grain surface regresses according to 
Figure 4-1, where the inside corner at P produces the circular burning region 
pictured, due to normal regressions. This region appears after a distance 
equal to the length of the inhibitor buried m the grain has been burned. 



Figure 4-1. GRAIN REGRESSION IN A SLOT 


The slope of the gram at the slot face (foreface of the slot aft face of 
the propellant) can be measured from the values input to the reference planes. 
The surface can then be placed m the linear form 

y = xtan(f) + r^ - Atan^ + s^/cos<j> 

where A is the longitudinal position, as indicated, along the motor axis; ^ 

is the port radius at A; and s, is the distance burned. 

b 
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After burning has progressed so that the surface impinges directly on the 
circle, the following equations can be solved for the intersection point: 

2 2-2 
circle: (-x - A-) + (y - f )' = (s^ - T) 

line : y = xtan$ + r - Atan<j> + s^/cos<|> 

where r^ is the distance from the centerline to P, at A, and T is the thickness 
of the inhibitor buried m the grain (Figure 4-2) . 


If the line equation is solved for s fa/ the distance burned, then a form 
appears such as 

s^ * cos<j>(y - r^ + Atan<}>) - xsin«jj 

Taking the coefficients of x and y, and noticing that 

2 2 
cos <f> + sin <j> = 1 

it can be stated that the intersection tracks out a parabola, again called a 
bum-back parabola. 


Substituting the value of s^, in terms of x and y into the circle equa- 
tion, obtain a form such as 

Ax 2 + Bxy + Cy 2 + Dx + Ey + F = 0 B 2 - 4AC = 0 


where 


A = cos $ 


D = - 2 (A - sin<f>cos4> (r^ - Atan4> + T/costb)) 


B = 2sin$cosij> = sin2<ji E = - 2(r^ - cos <j> (r^ - Atan<£ +• T/cos<j>)) 


. 2 

C = sm $ 


2 2 2 2 
F = A + r 2 - cos <j> (r^ - Ataniji + T/cos<f>) 


where T is the thickness of the inhibitor imbedded in the grain. (Note that 
B 2 - 4AC = 0.) ' 


Using the quadratic formula, the above second^degree equation for x is 
solved in terms of y. Knowing this, and knowing the case radius, the point of 
intersection of the bum-back parabola and the case can be computed. This 
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point then becomes an important position, because the region from this point 

to the slot interface will be considered one mass addition region. Notice 

that this length has a limiting value, that being the length r - r where r 

C X c 

is the case radius at A (Figure 4-2) . 



Figure 4-2. AN ANALYTICAL VIEW OF A SLOT REGION 
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Since the burn-back parabola has an equivalent form m terms of x and 

instead of x and y, and since is known at any time, the position point R 

can be calculated at any s, . Port volume and burning surface area m this 

b 

region are calculated as volumes and areas of revolution. The formulas for 
the "arced" region are developed in Appendix B, while the usual cone frustum 
formulas are used for the linear region. In order to keep equation consistency, 
it is assumed that the gram regression is constant in this particular mass 
addition region, that is, the burn rate at the slot face is the same as it is 
at the point of intersection of the burn-back parabola with the case. This is, 
however, a very small price paid for analyticity in this region. 

At the beginning of the bum, the point R does not travel along a burn- 
back parabola, but along an angle bisector. In this region, special care has 
to be observed to include m the surface area calculations the annular region 
at the end of the grain (Figure 4-3) . That is, when region 1 is revolved about 
the motor axis, the annulus created has a definite burning area. It is calcu- 
lated as the difference m the areas of two circles, and added to the rest m 
that region. 



Figure 4-3. TRACKING THE PATH OF ABRUPT GEOMETRY CHANGE IN A SLOT REGION 
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The path that point R takes during the burn is also indicated in Figure 
4-3. Because of normal regression, R moves along the angle bisector of <)A^BP 
until the burning surface has reached the level of point P. Point R then 
travels along the burn-back parabola to the case. There is a smooth transition > 
from one curve to another, since the angle bisector is tangent to the burn- 
back parabola at Q. In effect, then, the path of R has a continuous derivative. 

In Figure 4-4, the burning surface has regressed to a point where the cir- 
cular region cuts the case insulation. Knowing the radius of the circle as 
s^ - T, and r^ at A, then the y-coordinate of point Q can be ascertained. 

Knowing, from input data, the slope of the grain at Q, the x-coordinate of $ 
is known. Since points P and S are known, angle <)OPS can be calculated, a 
further decrement of the circular sector established, and the true length of 
arc RS kept up with. Points Q and R then become the points (x^, y ) and 
(x^» y^) of Appendix B so that volume and burning surface area can be calcu- 
lated. 



Figure 4-4. USE OF A BURN-BACK PARABOLA IN A SLOT 
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Section V 

TRANSITION REGION 

The region of transition between the star section and the CP section in 
the case where there is no inhibitor between them presents a realm of diffi- 
culties hitherto unencountered. What burning exists in the rest of the motor 
is one dimensional (radial) at best and easily two-dimensional at worst. Exam- 
ples of one-dimensional burning exist in the CP section, and two-dimensional 
burning m the star and slot regions. The transition region, however, presents 
a third spatial dimension to the burning. 

The region now under consideration, here called the transition region, 
is that region, pictured m Figure 5-1, where complex geometry appears during 
the burning because there is an uninhibited CP face presented m each star 
valley (Figure 5-2) . 


<L~ 



Figure 5-1. THE LIMITS OF A TRANSITION REGION 
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Figure 5-2. INITIAL RELATIONSHIPS AT A STAR-CP INTERFACE 


The physical situation then exists where, during the burning, a newly 
created star region emerges from the CP region. 

The angle ct^ of Figure 5-1 is generally not a zero angle in large motors . 
A tilt in the direction indicated will facilitate mandril withdrawal, and so 
will usually be a positive, but small, angle measure. 

Figure 5-2 represents a view looking into a star valley in the transition 
region. Notice that it is the burning on face ABCD that produces the new star 
geometry. The radius of arc CE is generally not the fillet radius m the star 
trough. The object AFGHJIDA is due to the fillet radius and, produced 
longitudinally, closely represents a quarter of a cylinder. The region 
between D and I, however, since the cylinder is wrapped against the valley side 
on the one hand and the valley floor on the other, produces a toroid. Because 
of the slant produced by positive angle a the region AFGD is again a cylinder, 
but now slanted off the perpendicular. 
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At ignition, the surfaces represented in Figure 5-2 all move m a direc- 
tion perpendicular to the surface tangent at each point. As mentioned pre- 
viously, area ABCD moves longitudinally, creating a new star-shaped region. 

The plane of the star valley wall moves in its usual direction. The toroidal 
channel grows m radius. 

At a time further in the burn, a shape like that of Figure 5-3 emerges. 
Since the governing force that controls the star region (not the transition 
region) is the moving together of the two planes representing the valley sides, 
it can be seen that under the supposition of constant bum rate there will be 
a time after which the line joining the star peaks will be parallel to the 
motor axis. That supposition is not made; however, this line should vary 
from the parallel by only a small amount. After the line reaches parallelism 
(under the assumption, again, of constant burn rate) , the whole star section 
is of homogeneous geometry; that is, the geometrical mensuration quantities 
at any plane perpendicular to the motor centerline are the same. 



Figure 5-3. A STAR-CP INTERFACE DURING BURNING 
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While the general case is where the burn rate varies longitudinally 
along the motor axis, it is helpful and not too incorrect to assume constant 
bum rate. The computer code developed does not assume this, but the events 
themselves are intuitively easier to picture with that qualifier. 

Notice in Figure 5-3 that a new star region is still being created as a 
result of the regression of plane ABCD. This continues throughout the burn, 
the only change being the subtle change occurring when ABCD bums out. The 
toroidal channel of Figure 5-2 has grown m dimensions; its center, however, 
remains in the same position. In Figure 5-3, web bum-through has occurred 
in the star valley, as seen from the presence of a bare case. It will be 
noted that the star peak region in which point A' is placed is receding at an 
accelerated rate, dependent only on the case sector angle and the angle of 
the valley side. Care must be taken in any analytical routine to include this 
fact. 


Attention is now turned to Figure 5-4, which represents a time still 
further into the burning. Plane ABCD of Figure 5-2 is gone, as is most of 
the original star region. Notice that a star is still being created because 
of a longitudinal burning component at point B (Figure 5-4). That_burnmg is, 
of course, governed by the radius of the arc connecting the star valley floor 
and the forward face of the CP section. This radius grows at the rate given 
by the bum-rate, and proceeds until no propellant is left. The growth of the 
toroidal channel produces the curved valley side. It is also the driver for 
telling where complete grain annihilation has occurred. 

The modeling of the transition region is achieved by tracking seven 
important points, and then using a star geometry program based on the one 
described above, augmented to be more general. An increment-dividing plane 
is placed at each of these seven points; and as these points move, the IDP's 
associated with each point are allowed to move with them. Before going into 
specifics concerning the placement of these points and the philosophy asso- 
ciated therewith, time will be taken to describe the generalized star geometry 
model. 
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Figure 5-4. A STAR-CP INTERFACE LATE IN THE BURN 


As just discussed, the star geometry routine needs the following pieces 
of information {Figure 5-5) : 



initial port radius 
distance burned 

distance from centerline to bottom of star valley - r^ 

radius of the fillet m the star valley 

width of the star valley at the top 

width of the star valley at the bottom 

radius of the case (centerline to insulation) 


N . number of star points in a cross section 
stars 

All other necessary constants can be calculated from this set. In terms of 
how the original star routine used this data, it was noticed that the set is 
redundant. 
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Figure 5-5. QUANTITY DEFINITIONS IN A STAR SECTOR 


The terms involving r^, s^, and r f can be replaced with the two new terms 


r l = r i + S b 


r f = r f + S b 


that is, at any given time m the bum, the port radius is given by r^. But 

then, the values of r^ and need not be stored; and instead of obtaining 

values from r^ and s^, r^ can be used effectively with s^ set at zero. All 

information necessary to the routine can then be input as r, or r. + s . 

. _ 11b 

Similar statements cam be made concerning and s^. 


To illustrate this point, consider Figure 5-6. Since the geometry 

routines, in effect, work with present port radius and present fillet radius 

(this is a simplified approach, but accurate), then it doesn't matter whether 

the routine is fed r. and s, , or r„ and s, and calculates r„ and r" , or whether 
__ ibrb If 

r^ and r^ are fed directly as initial port radii and initial fillet radius 
and the bum distance is set equal to zero. The two approaches are the same, 
both logically and computationally. 
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Figure 5-6. NOMENCLATURE FOR AN UPDATED STAR MODEL 


Consider now the transition region ]ust described. There are some 
semantics necessary before proceeding to the seven moving points. In the 
star geometry routine, the value r^, either as input by itself or as input 
in the 'form r^, is the initial port radius. In the region under discussion, 
then, stands for that value of the port radius existing at a point the 
instant that the geometry at that point changed from CP configuration to a 
star configuration. Although this value is essentially not the initial port 
radius, it is the initial radius as concerns geometry calculations at that 
point as a star configuration. This sequence of events is shown in Figure 
5-7. At ignition, the situation in (a) prevails, with an initial port 
radius, r^. In (c) , the moving plane has progressed to the point where the 
star is being burned into the cross sectional plane at A. The time that 
occurs is represented in (b) . Notice that the advancing plane is directly at 
position A. At any time after this, the existing geometry at A is star shaped; 


at any time before, it is CP. That value of s^ which places the advancing 
plane right “at A is called s * , and the value + s* is then taken to be the 
initial port radius at A (at that time when the geometry at A changed radically). 
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Figure 5-7. FORMATION OF A STAR VALLEY IN A PREVIOUSLY CP SECTION 
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The subtleties of those distinctions are necessary to properly supply 
the generalized star geometry routine with the necessary input parameters . 

The specific geometry in this region points out why it is necessary to have 
the CP routine working on a sector basis — in order to have the necessary 
continuity between star and CP. The calculations for r can also be. made for 

r f 

5.1 TRACKING SEVEN POINTS 

In order to adequately model the transition region, the region must be 
broken into at least six mass addition regions. At the beginning of the burn, 
some of these regions are degenerate, while at the end, other regions are 
degenerate. (A degenerate mass addition region is one between two idp's 
occupying the same position. ) The seven positions bounding these six mass 
addition regions will be the points tracked to determine the transition region 
geometry. The seven points have been chosen either to simplify volume and 
surface area calculations or to indicate abrupt changes in port geometry. 

Figure 5-8 presents the initial position of the seven points , named 
to P_,, along with the path that each will take during the burn. The rest of 
subsection 5.1 deals with the definition of these points and the mathematical 
equations that govern their motions. 

5.1.1 Tracking Point 

Point P^ is established initially at the point just before the large taper 
which takes the star into the CP region. It marks this position until it 
reaches a point just before it moves down the arc indicated in Figure 5-8. 

When it reaches the arc, P^ moves down a burn-back parabola. While on this 
arc, marks a position of radical geometry change. 
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Figure 5-8. THE PATHS OF SEVEN IMPORTANT POINTS 


The reason for P-^'s strange descent before the arc is because that once 
the concave portion of a star peak is burned out, the star recedes at an 
increased rate. Figure 5-9 presents a more intimate view of the path of p 
From surface AO^ to surface BO^ (AO^ is the initial unburned surface) the 
point P^ progresses down an angle bisector. Since the port radius is input 
at the initial position of P^ through a reference plane there, and this same 
item is input for the two adjacent reference planes surrounding P^, then the 
slope of the grain between A' and C^, and between and A can be determined. 
Using the arc tangent function, angle a can be determined and so the angle 
°2°1 A iS known * P i proceeds along this angle bisector until the concave 
region of the star point at 0^ is burned out. 

When point P^ reaches 0 , the time begins where the grain regression is 
accelerated in region B'O^ However, the acceleration is just a constant 
multiple of the bum-rate, as determined by the procedure described next. 

I 

Suppose the burn distance 0 3 0 is given to be 1. It is desired to ob- 
tain the grain regression constant, R, which will be the distance 0 0^. The only 
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Figure 5-9. AN ANALYTICAL VIEW OF THE PATH OF P-, 


things that will govern R are the angle of inclination of the star valley side 
(angle <f) and the sector angle <P , that is, \ji = ir/n where n is the number of 
star points. By the choice of coordinate axes (Figure 5-2) the length R is a . 
length on the star peak centerline. 

The two line formulas are given by 
- Y - xtan$ 

y = - xtantj) + s^/cos^ + C 
where C is a constant 

Substituting 0 and 1 for s fc (the points represented by 0 2 and o' 2 
respectively) obtain 

y^ = - x^taniji + C 
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= - x 2 tan<J> + l/cos<f> + C, 


so that y 2 - y x = - (x 2 - x x ) tan<j) + cos<£ 


From the star peak centerline equation, obtain 


y 2 - y-L = (x 2 - x x ) tanij; 


so that 


(tan^ + tan(j>) (x 2 - x^) = l/cos<j) 


or 


X„ “ X, = 


2 1 cos<|> (tanijj + tan<J>) 


Translating this distance to the star peak centerline there is 

1 


R - 0 0' = 


x 2 - x x 


2 2 cos$ cos^coS-fi (tani|> + tanifi) 


As noted, R will be called the grain regression constant. 


The path that point P^ takes immediately after passing point 0 can now 
be described. For an analytical view of the situation, consider Figure 5-10. 
Side AB is the initial grain surface, and CE is the burning surface after 
burning length CB from AB. CB is set to 1 in the above example. Since angle 
a is known; by congruence methods, angle £GBF also has angle measure a. 

Knowing that £GBH is right, then £FBH =? tt/ 2 - a. From this, and opposite 

IT 

angles, £AB1 = ir/2 - a. Likewise, £BCH = — - a. From the above discussion, 
it is known that length AD = R, the grain regression constant, where BC = 1, 
so BH = cosa, CH = sino. But then, IB = R - cosa so that AB = (R - coscO/s-ina 


Hence , 


or 


tan t, = 


1 

R - cosa 
sma 


tan Xy = 


sma 


R - cosa 


But then, 

C 


arctan (— 
R 


sing 
- cosa 
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Figure 5-10. NOMENCLATURE TO DEFINE THE PATH OF P ] 

Knowing £, then angle £DA C is known, and the path of after reaching 
(of Figure 5-9) is known. It is stressed that point cannot be truly 
extended down the angle bisector of -jA'o^A because of this accelerated nature 
of the regression of the surface B'C^. In fact, the true path of is a 
radical divergence from the angle bisector, in general being in the opposite 
direction. 

is now partially tracked. For the remaining portion of its path, 
consider the situation of Figure 5-11. The surface of the grain is given by 
curve DBCE, the arc BC arising from the inside angle at A, and AB - AC. 

Point P now is at B. Notice there is a sharp peak at B, indicating an 
abrupt change in geometry. This situation does not exist at C since line CE 
is tangent, to arc BC at C. is then given the job of tracking the point B 
through time to be able to demarcate the change in geometry. At any time, 
arc BC is given by a circle equation and line DB given by a trivial line 
equation. The intersection of the two, through time, gives rise to a burn- 
back parabola, along which P^ will travel. The mathematics of the derivation 
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of the path of P^ in this section of the bum is like that done above and is 
omitted. 


p 1 



Figure 5-11. POINT P ] TRAVELS ALONG A BURN-BACK PARABOLA 


The mensuration quantities obtained for the mass addition region ending 
at P^ are obtained m the same manner as those m the star region. 

5.1.2 Point ?2 (Stationary) 

Point P 2 is included as an increment dividing plane location only, and 
exists so that the mass addition region between P^ and P 2 does not become too 
large. It is positioned at the center of the circle connecting the star valley 
floor with the CP face (Figure 5-12) . It is not allowed to move, and the 
geometry there is always star geometry until burnout occurs. 

The values of grain surface area and port volume between P^ and P 2 belie 
the simplicity of the movement of P^ . Because of the circular trough appearing 
at A in Figure 5-11, the grain surface is not linear. For this reason, and for 
the reason that there does not exist a surface area formula as general as the 
prismatoid volume formula, it was deemed necessary to cut the interval between 
P 1 and P 2 into several subintervals, calculate initial grain constants at each 
cut point using the star geometry routine, and calculate the volume and surface 
area of the mass addition region as the sum of those quantities of each sub- 
interval . The method used is to once again depend on the cone frustum formulas , 
which essentially means a linearization, or better, a piecewise linearization 
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of the grain between P 1 and The number of subintervals is given as 'a 

trade-off between accuracy and expediency. 



Figure 5-12. THE STATIONARY POSITION OF POINT ? z 


5.1.3 Tracking Point Pg 

P^ is used to mark the position of point A in Figure 5-4 — the point of 
grain annihilation. The point, of course, remains stationary until a given 
distance burned, and is then tracked in the aft direction. 

The governing geometric figures which give rise to the analytic equations 
governing the motion of P^ are the toroidal channel in the star valley and the 
plane which represents the centerline of a star peak produced longitudinally. 
As the toroid expands during the burning, it will intersect the mentioned 
plane. At the time of first intersection, the entire valley wall of the star 
has been eliminated. When the toroid gets large enough to reach the line 
given by the plane's intersection with the case wall, total annihilation of 
the grain at that point is assured, once the toroid reaches this line, the 
point of intersection recedes toward the aft end. This point is P^. 

Using the Cartesian coordinate system described m Figure 5-2, the 
equation of the toroid (with respect to subsection 2.2) is given as 

(x - r 2 ) 2 + (z - P) 2 = + s b^ 2 ~ <Y - r 4 ^ 2 + r 2 " r 2 ) 2 

where 
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and are described above, and 
P = longitudinal position of plane P 2 



= radius of the circle connecting the valley floor with the 
CP face 

r„ = r, - r_/cos$ 

4 hot f 

<j) = angle the star side makes with the x-axis , 

The line given by the intersection of the centerline plane with the case wall 
xs given in three space as 

x = r cosi p 
c 

y = r c sinif; 
z 

Substituting this line into the toroid gives 

V | 2 2 2 2 

( f (r. + s, ) - (r siniji - r ) + r, - rJ - (r cos!{i - r ) 

1 f b c 4 of c 2 

where P is the position (longitudinal) of P and s. > s. ,where s, is the 

2 b — b ,m b ,m 

maximum bum distance at position P 2 - Everything in this equation is known, 
and so can be determined directly; that is, the point of total grain 
annihilation can be calculated at any time during the bum. 

The point is initially positioned at position P 2 , and is stationary until 

s, > s, . It is then allowed to progress toward the aft end (determined by 
b — ■ b,m 

the choice of the plus sign in the above expression giving z) . The geometric 
mensuration quantities existing at P^ are simply made by the following 
instructions. ' 

(1) If s, < s. quantities at P, = quantities at P„ 
b — b,m 3 2 

1 2 

(2) If s, > s. port area is maximum at — r i p, port perimeter 
b — b,m 2 c * 

is 0. 
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When s j-, s b m , then there is a positive distance separating P 2 and P 3 - 

Since no grain exists between these two planes, by the very definition of P^ , 

then no grain surface exists, and the port volume calculations are made using 

the frustum formula described above. To induce generality, notice that this 

same situation holds when s, < s, , since then the distance between P„ and P_ 

h ~ b,m 2 3 

is zero, and so no surface area exists and the frustum formula for volume can 
be used, but will also yield zero. 

5.1.4 Tracking Point P 4 

P^ marks, again, a position of abrupt geometry change. It is also 
placed at position P 2 at the beginning of the burn, and is used to mark the 
burnout of the star valley floor (see the region about point F in Figure 5-3) . 
Its position is controlled by the growth of the circle (cylinder, actually), 
whose center is marked as P 2 m Figure 5-12. Point remains stationary at 
position P until s > s, where s. gives the web thickness in the valley 

2 JO ~~ D/V P / V 

floor. Once this circle reaches the case wall, the growth of the circle then 
forces P^ toward the end of the segment. 

In the Cartesian coordinate system described in Figure 5-2 the circle and 
case plane are given by 

circle: (x - r 2 ) 2 + (z - P) 2 = (r 3 + s b ) 2 

plane: x = r^ 

where P represents the longitudinal position of plane P 2 - Upon substitution, 
obtain 

(r c - r 2 ) 2 + (z - P) 2 = (r 3 + s b ) 2 
or 

where the plus sign denotes movement toward the aft end. 

When s b <_ s b - ,the cross sectional port area and port fuel perimeter are 
assigned those at P 2 . when :> s fa v ,those values are obtained through the 
use of the modified star geometry routine. The volume and grain surface area 
between P 3 and F ^ are obtained by cutting the distance between P 3 and P^ into 
halves, using the cone frustum surface area formula m each half, and using 
the general prismatoid volume formula over the whole region. 
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5,1.5 Tracking Point Pg 

The next point to be tracked is the point which specifies the position of 
burnout of the concave portion of the star peak. This is position Pg on 
Figures 5-4 and 5-8. From F-igure 5-2, it is clear that the quarter cylinder 
AFGD will control the position of P<. for some of the burn, and the torus at 
the bottom of the valley will control it for the rest. 

At the beginning, Pg is placed at position A of Figure 5-11. It remains 
at that point until the concave portion of the star peak has burned out. In 
effect, Pg is placed at A until the cylinder has grown large enough to inter- 
sect the star peak center plane. It is then allowed to travel in the positive 
z direction (Figure 5-2 for coordinate reference) . When the burning surface at 
P^ has reached a level close to then the motion of Pg is governed by the 
toroid mentioned above. This condition exists until Pg reaches the case wall. 

Notice from Figure 5-1 that the cylinder leans m the positive z direction 
by angle a^. It also leans in the positive y direction by angle <f>, the slant 
of the star valley side. In accordance with subsection 2.2, the cylinder has 
equation 

[ 

If ctg is the slope of the grain immediately after the initial position of 

Pg, then the port itself is a cone with equation 

2 2 2 
x + y - (tanctgZ + r*) = 0 

where r* is the port radius at the initial position of Pg . At any particular 
time during the burn, the port surface, itself, has equation 

x 2 + y 2 - (tanagZ + r* + s^) 2 = 0 

where is understood again, to be the burn distance. The annihilation takes 

place at the star peak center plane which has equation 

x = y/tanifi. 

Substituting this into the cone equation, obtain 


y - (r 4 + z cos f— 
r f + s b 


z sinC^” + <j>) cota^) 


(r f + s b )/sina 4 
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y = (tanc^z + r* + s^) simp 


Hence, 

x ~ (tana^z + r* + s^) cosip 

Using these values and substituting back into the cylinder equation, there is 

— 3tt 

(tana,z + r* + s,) cosip - (r_ + sinOr— + <f>) cota.z) 

3 J. D 2 2 4 

(r f + s b )/sxna 4 

(tana^z + r* + s b )sinip - (r^ + cos + <f>) cota^z) 

r f + s b 

which will reduce to a farm such as 
2 

az + bz + c = 0 




where 

2 2 2 
a = {tana^simp - sinipcota^} + sin (tana^cosip + cosipcota^) 

b = 2 [ (tana^s imp - simpcota^) (sinip (r* + s^) - r^) + 

2 — 
sin (tana^costp + cos^cota^) (cosip (r* + s^) - r^U 

2 2 —2 2 
c = (simp (r* + s ) - r ) + sm a (cosip (r* + s ) - r.) - (s + r ) 

1 b 4 4 lb2 bf 

and so the position of P^ can be deduced using the quadratic formula. 


After the burning has taken the surface to about the level of , then 
the motion of P 5 is governed by the intersection of the torus with the star 
peak center plane. 


As seen, the torus has equation 

(x - r 2 ) 2 + (z - P) 2 = (y(r f + s b ) 2 - (y - r 4 ) 2 + r 3 - r f ) 2 

The plane has equation 

x = (ztanct„ + r* + s„ ) cosip 
3 1b 

y = (ztana_ + r* + s,) sinip 
3 I b 

Substituting these expressions for x and y above, yields 
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— 2 2 

(ztanc^cosip + (r£ + s^) cos \p - r^) + z = 

(\{r f + s b ) 2 - (ztana 3 sini(i + (r* + s b > sintj; - r 4 ) + r 3 ” r f ) 2 

Notice that this is a fourth degree equation in z. Because of the compli- 
cated nature of this particular equation, a Newton-Raphson iteration process is 

utilized, with an initial guess of z n - s , the process converges rapidly. It 

o b 

is estimated that in this case only four iterations are necessary to bring the 
answer to within seven decimal place accuracy. 

The mensuration quantities at P_ are obtained by use of the modified star 

D 

geometry routine. The surface area and port volume between P . and P are 

4 5 

calculated in the same way as those between P„ and P„. 

3 4 

5.1.6 Tracking Point Pg 

The path of P c is probably the simplest, other than that of P„. Notice 
o 2 

that Pg begins beneath point A (Figures 5-8 and 5-11) at the point of tangency 
of the CP face with the circle in the star valley. It progresses in a linear 
fashion until it overtakes point at which time the two are considered the 
same. 

The governing equation is 

z = s, cosa. 
b 4 

x = r. + r,sina. + s^sina„. 

2 3 4 b 4 

The mensuration quantities here are obtained from the modified star 

geometry routine, and the port volume and surface area are obtained in the same 

manner as that between P_ and P . . 

3 4 

5.1.7 Tracking Point Pj 

The final point to be tracked is point P^. At any time, points P ? and P 1 
mark the boundary of the defined transition region. Point P_, is initially 
placed at the intersection of the CP region with the star at the aft position 
of the taper of a CP face (Figure 5-8). At the beginning, P ? is tracked along 
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the angle bisector as shown in Figure 5-8. The angles are all known, so the 
angle bisector is known. The governing equations are 


. . a. + a. 

s in/ 1 4 3 


z* + s. 


\ 4 


sm { 3 4 ) 

V 4 2 ' 


x = (z - z*> tanct^ + s^/cosa^ + x* 


where z* is the initial position of P ? and x* is the port radius, there. 


When Point P- catches P_, as it must since they travel at the same rate 

O / 

but move in different angles, then point P^ must travel down the burn-back 
parabola evolving from the grain surface and the intersection of the circle 
growing from position P 2< that is, the burn-back parabola generated from the 
equation set 

/ . — 2 2 2 
circle: (x - r_) + (z - P) = (s, + r_) 

1 Do 

line: x *> (z - z*)tana_ + s, /cosa_ + x* 

3 b 3 

where P is the longitudinal position of point 


The mensuration quantities are calculated by the CP routine, and the port 
volume and surface area of the region between P^ and P ? (when they are not 
zero) are calculated m the same manner as that between P^ and P 4 - 

5.2 SUMMARY OF TRANSITION REGION 

In Section V, the quantities of port volume and port fuel surface area 
were calculated as corollaries of the positions of seven points . These 
formulas ultimately depend only upon burn distance, and the particular depen- 
dent variable is calculated anew at each new value of burn-distance. In this 
manner, any computational errors are not accumulative, but should remain at the 
same level. 


The difficulty of the region arises from the appearance of nonlinear star 
sides. This forces linear approximations to be made because of the absence in 
the literature of generalized surface area formulas to handle nonlinearity. 
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Section VI 

GENERAL SUMMARY 


This document was designed to present an analysis of the geometrical 
changes occurring in a burning solid rocket motor. From very simple assump- 
tions about the grain regression, a unified functional approach was built to 
look at the regression problem as one of being functional in nature, depending 
only on the burn distance. 

The burning in each particular geometrical' region contributes to the 
overall result. The solution at any particular burn distance is initiated by 
specifying a bum rate and a delta time at each increment dividing plane. No 
assumptions were made as to the behavior of the burn rate, and it is allowed 
to vary in any fashion. 

Once the bum rate and delta time are established at each IDP , the code 
progresses from plane to plane, calculating port perimeter and port area. From 
these quantities, the numerical values of port surface area and port volume are 
calculated for the mass addition region between the two planes. In a slot, ' 
only port volume is needed since no fuel surface area exists. This procedure 
is continued until the aft tangent plane is reached, at which time the code is 
ready for the next case. 

Information from the reference plane input is used by the code to deter- 
mine where the geometry delineations occur. 

The results of this task are twofold: firstly, the incorporation of this 

geometry model into the SRB prediction program. No difficulties in this respect 
are foreseen; however, this incorporation is a large task in itself because of 
the tediousness and sensitivities involved. So that the reader can have a 
better understanding of what this code actually produces as useable quantities. 
Appendix C presents a sample case. Secondly, through equations like those 
presented here, and others depending on this geometry model, it is now felt that 
the theory exists to construct a computer program to calculate moments of inertia 
and centers of gravity of a burning SRM in a self sufficient manner. 
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APPENDIX A 

SUPPLEMENTARY GEOMETRY 
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THE BRETSCHNEIDER— HERON FORMULA* 


Consider the object m Figure A-l. 


A 



Figure A-l. A GENERAL QUADRILATERAL 


The Bretschnexder formula for the area, K, of quadrilateral ABCD is given 
by 

K = \j (s-a) (s-b) {s-c) (s-d) - abed cos 2 ? 

, A + C B + D , 1, 

where Z is either — - — or — - — and s = — (a + b + c + d), where it is assumed 

that the quadrilateral is convex (that is, no corners greater than 180 degrees, 

and not a skew projection from 3-space) . 

In order to substantiate the formula, consider the augmented form given 
in Figure A-2, where line segment BD has been drawn in and labelled e. 

Now, e itself is known by the law of cosines in two different manners: 

2 2 2 

e = a + d - 2ad cosA 


This appendix has been provided to correct an error given in the definition 
of the Bretschneider Formula in the 21st Edition of the CRC Mathematical Rand- 
book. This appendix provides the corrected equation and shows the theoretical 
development . 
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A 



Figure A-2. A TRIANGULARIZATION OF A GENERAL QUADRILATERAL 



and 
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Hence , 


2 2 2 
K — K 4- 7TC TC 4- 

□ ABCD AABD AABD ABCD ABCD 


1222 222 1 
= — [a d sm A 4- b c sxn C] 4 - — abed sinA smC. 


Using the cosine half Angle formula (remembering that no angle in the quadri- 
lateral is greater than 180 degrees) , and the formula for the cosine of the sum 
of two angles , obtain 


sinAsmC = 14- cosAcosC - 2cos 


2 A + C 


Then 


2 1 . 2„2 . 2 ,2 2 2 . 1 2 (A + C 

^□ABCD ~ "4 ^ Sln A + b c sin C] 4 - — abed [1 4 - cosAcosC - 2cos I — - — 


1 r 2 2 2 . ^ 2 2 
- [a d (1 - cos A) 4- b c (1 


- cos C) 4 - 2abcd 4 - 2abcd cosAcosC] 


- abed cos 


2 A + C 


2 2 

using sin a 4- cos a = 1. 


Factoring, obtain 


^dABCD ~ ^4 (^ a< ^ + hc)^ - (ad cosA - be cosC)^) - abed cos^ 


A 4- C 


Using the expressions for e above, 

12 2 2 

ad cosA » — (a 4- d - e ) 

12 2 2 

be cosC = — (b 4- c - e ) 

so that 

1 2 2 2 2 

ad cosA - be cosC = — (a 4-d - b - c). 


r/-vroTT TTY OF ® 
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Then, using this identity, it is easy to show that 

((ad + be ) 2 - (ad cosA - be cosC) 2 ) = (s-a) (s-b) (s-c) (s-d) . 


From this, the Bretschneider formula is obtained as 


But 


^bABCD ~ ^ S " b ^ (S_C ^ “ abed cos 2 ^ *" - ^ 


2 \ + C\ 2 (B + D\ _ . . . , 

cos | — - — I = cos I — 7 — , so the formula is stated 


k oabcd 


- \[{^ 


a) (s-b) (s-c) (s-d) - abed cos £ 


When one side of the quadrilateral is degenerate, that is, has aero length so 
that two vertices coincide, then another familiar formula arises. To see this, 
suppose that d » 0, that is, side AD is degenerate. 


Then using the Bretschneider formula 

K OABCD~ ^ s ~ Sl1 < s “b) (s-c) (s-0) - abc(O) cos 2 £ 
or 

Wd * \Z ,ta - a) |s_b) (s -°> 

which is Heron's formula for the area of the general triangle with s = ~ 

1 2 
(a + b + c + d) = 2 * (a + b + c) . 

Notice, however, that only Bretschneider* s formula need be referenced, 
for with a degenerate side, Heron's formula is obtained.* Referencing the star 
burning section, this property is used heavily, and the Bretschneider formula 
occurs with the name Bretschneider-Heron formula to denote the case when the 
formula is used to calculate the area of a triangle or a quadrilateral, or 
both. The two cases are shown simultaneously in Figure A-3. 
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Figure A-3. THE USE OF THE BRETSCHNEIDER-HERON FORMULA 


Knowing angles R, S, U, V, and sides RS, SU, UV, and VR, the Bretschneider- 
Heron formula can be used to calculate port area m quadrilateral RSUV. Further 
in the burn, R progresses to R' and sides R'S, SU, UR* are known; and the 
Bretschneider-Heron formula is used to calculate port area contributed by tri- 
angle R'SU, not mattering what the angles £UR'S, £r'SU, and £SUR' are. 

Attention now is turned to another star-burning anomaly, that being to 
calculate the area of an angular segment of a sector of a circle. Consider 
the situation described m Figure A-4. 

The objective is to calculate the area ABC. Angle a is known and 
0 = ir - a, so sin0 - sin(ir-a) = sma. 
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Figure A-4. THE DETERMINATION OF A PARTIAL SECTOR AREA OF A CIRCLE 
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Hence , 


the area of the pie-shaped piece ABC is given by 

* \ <a + b,2v - Wffi’ " 2 e ° Sln “ 

- i (a + b) 2 sin' 1 (esina/a) - ^ OEDOB , - J ec 


sina. 


This particular expression is a close approximation to the real case, 
since the term a qBD q B . is approximation to the true area of BDOB ' , which 
is not a true quadrilateral. 
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APPENDIX B 

SURFACE AREA AND VOLUME OF A TOROIDAL CHANNEL 
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P 


MENSURATION QUANTITIES IN A TOROIDAL CHANNEL 


Consider the situation of Figure B-l where two circular arc lengths are 
established in the x-y Cartesian coordinate system, governed from the circle 
formula 

(x-a) 2 + (y-b) 2 = r 2 . (B-l) 


X = a 



This formula can be replaced by an equivalent form by solving for y m 
terms of x as 


y 


+ \J 


r - (x-a)' 


(B-2a) 


which, of course, is not a functional relationship. Interest here is in the 
upper semicircular region, and a true functional relation appears as 


y 


b 


(x-a) 2 


(B-2b) 


From (B-l) , the derivative of y can be obtained implicitly as 

B-2 
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dy _ x-a 

dx ~ y-b 

so that 



(x-a) 2 _ (x-a) 2 + (y-b) 2 
(y-b) 2 (y-b) 2 


2 

r 


(Y-b 2 ) 



(B-3) 


(B-4) 


The length of arc can now be calculated for the arc between and x^ (on 
either side of the line x = a) as 




If the arc from x_ to x (m either case) is revolved about the x-axis, a 

^ X 

toroidal channel emerges. The surface area and volume of the channel will 
now be calculated (Figure B-2) . 


The general formula for the area of revolution is, or course. 
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Figure B-2. A TOROIDAL SURFACE OF REVOLUTION 


which interprets, in this case, as 



using (B-2b) and (B-4) . 
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= 2irr (x 1 - x ) + 2rrbr 


r *1- 

arcsxn 

L r 


arcs in 




• arc length - 


= 2ir(r(x^ - x^) + b (arc length)). 


The general formula for finding the volume inside a surface formed by 
revolving a planar arc length about the x-axis is given by 


V — it I Tfy dx. 


Adapting this general formula to the particular case at hand yields 


x i~ — i. 

I* 1 i 2 2 * 

= tt j b+yr - (x-a) 


tt fb 2 + 2b \Ir 2 - (x-a) 2 + r 2 - <x-a) 2 |dx 


=* IT b 2 (x i - x 2 ) + 2b j 

x 

Xi _ 

- j (x 1 - a) 2 dx 


|b 2 (x x - x 


x^) + 2b \jr 2 - (x-a) 2 dx + r 2 (x 1 - x^) 


2 ) + r (x 1 - x 2 ) + b 


K- 


a) V r 2 - (x i - a) 2 


2 x^ - a f*2 2 " 2 X 2 ~ a 

+ r arcsm (x„ - a) \lr - (x„ - a) - r arcsm — — 

r 2 * 2 r 


(x - a) 3 (x_ - a) 3 

± H — — 

3 3 
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APPENDIX C 
SAMPLE RESULTS 
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